Strain controll of magnetic anisotropy in (Ga,Mn) As microbars 
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We present an experimental and theoretical study of magneto crystalline anisotropics in arrays of 
bars patterned lithographically into (Ga,Mn)As epilayers grown under compressive lattice strain. 
Structural properties of the (Ga,Mn)As microbars are investigated by high-resolution X-ray diffrac- 
tion measurements. The experimental data, showing strong strain relaxation effects, are in good 
agreement with finite element simulations. SQUID magnetization measurements are performed to 
study the control of magnetic anisotropy in (Ga,Mn)As by the lithographically induced strain re- 
laxation of the microbars. Microscopic theoretical modeling of the anisotropy is performed based 
on the mean-field kinetic-exchange model of the ferromagnetic spin-orbit coupled band structure of 
(Ga,Mn)As. Based on the overall agreement between experimental data and theoretical modelling 
we conclude that the micropatterning induced anisotropics are of the magnetocrystalline, spin-orbit 
coupling origin. 
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I. INTRODUCTION 

Dilute moment ferromagnetic semiconductors, such as 
(Ga,Mn) As, are favorable systems for studying and utiliz- 
ing controllable magnetic anisotropy since micromagnetic 
parameters of this ferromagnet are very sensitive to Mn 
doping, hole concentration, lattice strains, and temper- 
ature. The magnetic moment density is small in these 
ferromagnets and therefore the spin-orbit coupling in- 
duced magnetocrystalline anisotropy typically dominates 
the dipolar-field shape anisotropy. 

The control of the magnetocrystalline anisotropy in 
(Ga,Mn)As epilayers has been achieved by choosing dif- 
ferent substrates and therefore different growth induced 
strain in the magnetic layer, by varying the growth pa- 
rameters of the (Ga,Mn)As film, and by postgrowth 
annealing. 1 2 Reversible electrical control of the mag- 
netocrystalline anisotropy has been demonstrated by 
utilizing piezo-electric stressor!^ or by electrostatic 
gating in thin-film (Ga,Mn)As field effect transistor 
structures P^l Recently, a local control of the magne- 
tocrystalline anisotropy has been reported, which pro- 
vides the possibility for realizing non-uniform magneti- 
zation profiles and which can be utilized, e.g., in studies 
of current induced magnetization dynamics phenomena 
or non- volatile memory devices. 8 9 In these studies an ef- 
ficient method of local strain control has been used which 
is based on lithographic patterning that allows for the re- 
laxation of the lattice mismatch between the (Ga,Mn) As 
epilayer and the GaAs substrate.™^ The modification 
to the strain distribution can cause strong changes of the 
magnetic anisotropy for strains as small as 10 -4 . The 
high efficiency and practical utility of the lithographic 
pattering control of magnetic anisotropy in (Ga,Mn)As, 
demonstrated in the previous works, have motivated our 
thorough investigation of the phenomenon which is pre- 
sented in this paper. Our study is based on combined 



high-resolution X-ray diffraction and magnetization mea- 
surements and on macroscopic modeling of the strain re- 
laxation and microscopic calculations of the correspond- 
ing magnetic anisotropics. 

We investigate two sets of lithographically patterned 
(Ga,Mn)As microbars which differ in the thickness to 
width ratio, Mn doping, and hole concentration. First, 
we study the structural properties by high resolution 
X-ray diffraction of microbars patterned in the thicker, 
higher Mn doped as-grown (Ga,Mn)As material which 
has a large growth induced strain. The spatial distribu- 
tion of the lattice relaxation in the stripe cross-section is 
determined by comparing the measured intensity maps to 
maps simulated using the theory of elastic deformations 
and the kinematic scattering theory. The good agree- 
ment of the measurement and simulation shows that the 
applied model is quantitatively reliable in predicting the 
local lattice relaxation in patterned epilayers subject to 
small lattice mismatch. This allows us to infer the much 
weaker lattice relaxation in stripes fabricated in the thin- 
ner and lower Mn concentration (Ga,Mn)As by perform- 
ing only the elastic theory simulations. 

In the next step, we measure the magnetic properties of 
our samples by Superconducting Quantum Interference 
Device (SQUID) and extract the anisotropy coefficients. 
Stronger focus is on stripes fabricated in the thinner, 
annealed (Ga,Mn)As epilayer where the SQUID magne- 
tometry data allow for a reliable extraction of the tem- 
perature dependence of the anisotropy coefficients and 
for direct comparison with the microscopic model. We 
assumed a linear superposition of the in-plane uniaxial 
anisotropics and the presence of a single magnetic do- 
main when analyzing the SQUID magnetometry data. 
We show that the easy axis can be rotated by 90° by the 
micropatterning, completely over-writing the underlying 
material anisotropy at all studied temperatures. 

Finally, we calculate the anisotropy coefficients for 



a range of material parameters and temperatures be- 
low Tc> The lattice relaxations determined form the 
X-ray diffraction measurement and from finite element 
simulations are the inputs of the microscopic calcula- 
tions of the magnetocrystalline anisotropy. The micro- 
scopic model we use is based on an envelope function 
description of the valence-band holes and a spin repre- 
sentation for their kinetic-exchange interaction with lo- 
calized moments on Mn 2+ ions, treated in the mean-field 
approximation IMlMSl 



II. SAMPLES 

We study two sets of patterned (Ga,Mn)As epilayers 
grown on GaAs substrate. The samples in set A are 
doped nominally to 5% of Mn, annealed for approxi- 
mately 75 minutes at 180° C, and the epilayer is 25 nm 
thick. The Curie temperature Tq ~ 120 K corresponds to 
optimal annealing of the wafer. 16 The control sample Ao 
was not patterned. Samples Am^ and A[ 110 ] were pat- 
terned into 25 mm 2 arrays of stripes at an angle a ~ 140° 
and a « 50°, respectively. Here the angle a is measured 
from the [100] crystallographic direction. The uninten- 
tional 5° misalignment from the crystal diagonals caused 
by the microfabrication is accounted for when analyzing 
the data. The stripes are 750 nm wide, 100 jivn long, 
and separated by 450 nm gaps, as measured by Atomic 
Force Microscope (AFM). The fabrication was done by 
electron beam lithography and wet chemical etching us- 
ing a solution of phosphoric acid and hydrogen perox- 
ide. The AFM measurements revealed an etch depth 
of ~ 60 nm, and cross-sectional Scanning Electron Mi- 
croscope (SEM) imaging confirmed that the wet etching 
leads to anisotropic stripe cross-sections, with the A[ 110 ] 
stripes being undercut and the Amj i stripes overcut, as 
shown in Fig. [I] 

The samples in set B are doped nominally to 7%, not 
annealed, the epilayer is 200 nm thick, and the Curie 
temperature Tq ~ 85 K. The control sample B was not 
patterned. Samples Bm^ and B[ io] were patterned into 
arrays of 1 /im wide stripes with 1 /im wide gaps along the 
[110] and [010] crystallographic directions, respectively. 
The fabrication was done by electron beam lithography 
and dry chemical etching with an etch depth « 700 nm 
(B riioi ) an d ~ 900 nm (B[ io])- The sides of the stripes 
are slightly overcut in both cases owing to the symmetric 
dry etching. 

With respect to our theoretical modelling of the mag- 
netic anisotropics of our samples, we recall that relating 
the prediction to the measurement based on the material 
parameters is not straight forward due to the presence 
of unintentional compensating defects in the epilayers. 
Most importantly, a fraction of Mn is incorporated in in- 
terstitial positions. These impurities tend to form pairs 
with MnQa acceptors in as-grown systems with approxi- 
mately zero net moment of the pair, resulting in an effec- 




FIG. 1: (color online) Cross-sectional scanning electron mi- 
croscope (SEM) images of the stripes in set A. (a) Image of 
sample A[ 110 ] showing both the cleaved face and the top sur- 
face. Although difficult to discern, the profile is undercut. 
The curvature is due to the sample stage drifting during the 
exposure of the image. Introduction of coordinates fixed to 
the crystallographic axes and dashed coordinates fixed to the 
stripe geometry: the relaxation direction perpendicular to the 
stripes, the x axis, is rotated by angle uj — 45° with respect 
to the [100] crystallographic direction, the x axis. The angle 
uj = a — 45° describes the rotation of x with respect to the 
[110] axis, (b) Image of sample Ar^i showing a cut through 
the stripes and substrate in the x — z plane revealing the 
overcut sides of the stripes. 



tive local- moment concentration x e ff = x s — Xi. 16 Here 
x s and Xi are partial concentrations of substitutional and 
interstitial Mn, respectively. We emphasize that in fig- 
ures presenting calculated data the Mn concentration la- 
belled as x corresponds to the density of uncompensated 
local moments, i.e., to x e ff. 

Another input parameter of the theoretical modeling is 
the lattice mismatch which is different in set A and B as 
it depends on the partial concentrations of Mn atoms in 
substitutional and interstitial positions in the lattice and 
of other unintentional impurities.^ The lattice mismatch 
is determined by direct X-ray measurement as detailed 
in the following section. 

Fig. [I] introduces the coordinate system fixed to the 
crystallographic axes: x-axis along the [100] direction, y- 
axis along the [010] direction, and z-axis along the [001] 
direction which is the frame of reference for the micro- 
scopic magneto-crystalline anisotropics. The dashed co- 
ordinate system is fixed to the stripe geometry: x'-axis 
lies along the relaxation direction transverse to the stripe, 
y'-axis along the stripe, and z'-axis along the growth di- 
rection coinciding with the z-axis. The dashed coordi- 
nates are the natural reference for the macroscopic lattice 
relaxation simulations. 



III. LATTICE RELAXATION 

The lattice of thin (Ga,Mn)As films grown epitaxially 
on GaAs substrates is strained compressively due to a 
lattice mismatch eo = (a s — ao)/ao < where a s and 
clq are the lattice constant of the substrate and of the 
relaxed free-standing (Ga,Mn)As epilayer, respectively. 
The narrow stripes allow for anisotropic relaxation of the 
compressive strain present in the unpatterned epilayer. 
An expansion of the crystal lattice along the direction 
perpendicular to the bar occurs while the epilayer lattice 
constant along the bar remains unchanged. Parameters 
sufficient for determination of the induced strain are the 
lattice mismatch eo and the shape of the stripe, mainly 
the thickness to width ratio of the stripe. In the regime of 
small deformations the components of the induced strain 
are linearly proportional to the lattice mismatch. The 
strain tensor in the coordinate system fixed to the stripe 
reads: 
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where the lattice relaxation is quantified by p{x',z') 
which varies over the stripe cross-section, c\i and en 
are the elastic moduli. The strain components in this 
work are expressed with respect to a relaxed free-standing 
(Ga,Mn)As epilayer. Note that for p = the strain ten- 
sor e r describes the growth strain of the unpatterned epi- 
layer. In this section we investigate experimentally and 
theoretically the geometry of the stripes, the size of the 
lattice mismatch and the spatial dependence of the lat- 
tice relaxation p(x' ', z'). The results are used as an input 
of the microscopic modeling of the magnetic anisotropics 
in Sec. El 

Microbars in set B have larger thickness to width ra- 
tio than microbars in set A. Therefore the relaxation is 
expected to be larger in set B. At the same time, the 
(Ga,Mn)As epilayer has larger volume in set B, primarily 
due to a larger number of interstitial Mn in this higher 
doped unannealed material. The larger film thickness 
and larger growth strain in set B make these materials 
significantly more favorable for an accurate X-ray diffrac- 
tion analysis of the strain profile in the patterned micro- 
bars. 

The lattice relaxation in samples Br^ and B[ 110 ] was 
measured by high-resolution X-ray diffraction using the 
synchrotron source at ESRF Grenoble (beamline ID10B, 
photon energy 7.95 keV). For a reliable determination 
of both in-plane (v! x ) and vertical (u' z ) components of 
the elastic displacement field we measured the reciprocal- 
space distribution of the diffracted intensity around the 
symmetric 004 and asymmetric 404 reciprocal lattice 
points. The asymmetric diffraction was chosen so that 
the in-plane component of the corresponding reciprocal 
lattice vector h was perpendicular to the stripes. The 
diffracted radiation was measured by a linear X-ray de- 
tector lying in the scattering plane. 
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FIG. 2: (Color online) The measured (upper left panel) and 
simulated (upper right panel) reciprocal-space maps in the 
symmetric 004 diffraction of sample Bpio]- In the bottom 
row, the measured (points) and simulated (lines) intensities 
integrated along the horizontal (left) and vertical (right) di- 
rections are plotted. In the intensity maps, the color scale is 
logarithmic. 
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FIG. 3: (Color online) The measured (upper left panel) and 
simulated (upper right panel) reciprocal-space maps in the 
asymmetric 404 diffraction of sample B[ io]- In the bottom 
row, the measured (points) and simulated (lines) intensities 
integrated along the horizontal (left) and vertical (right) di- 
rections are plotted. In the intensity maps, the color scale is 
logarithmic. 



Figs. Inland [3] present examples of the measured (upper 
left panels) and simulated (upper right panels) recipro- 
cal space maps, showing two maxima corresponding to 
the reciprocal lattice points of the GaAs substrate and 
the (Ga,Mn)As layer. The bottom panels show the mea- 
sured and simulated integrated intensities for two direc- 
tions in the reciprocal space. Since the lateral stripe pe- 
riod was larger than the coherence width of the primary 
radiation, different stripes were irradiated incoherently, 



so that the lateral intensity satellites stemming from the 
lateral stripe periodicity could not be resolved. The mea- 
sured intensity distribution is therefore proportional to 
the intensity scattered from a single microbar. 

We fitted the measured intensity maps to numerical 
simulations based on the kinematic scattering theory and 
the theory of anisotropic elastic medium. We used a 
finite-element simulation (implemented in Structural Me- 
chanics Module of Comsol Multiphysics, standard partial 
differential equation solver) to obtain the local relaxation 
distribution p(x' ', z') in the stripes and derived the corre- 
sponding reciprocal space map. The angle of the sides of 
the stripes and the lattice mismatch eo of the (Ga,Mn)As 
and GaAs lattices were the two fitting parameters. The 
left column of Figs. [2] and [3] shows the measured diffrac- 
tion maps and projections. The right column shows the 
simulated results. The lateral and vertical projections 
of the measured and simulated intensity maps as well as 
the whole maps are used in the fitting. The coordinates 
q' x and q' z span the reciprocal space conjugate to the real 
space with coordinates x' and z' fixed to the stripe. They 
are measured with respect to the reciprocal lattice point 
004 and 404. 

The remarkable agreement of the measured and simu- 
lated diffraction maps shows that our model of the lattice 
deformations is quantitatively relevant in determining 
the local lattice relaxation p{x',z') in the stripes shown 
in Fig. |4j the lattice mismatch between the epilayer and 
the substrate, eo = —0.38 ± 0.03% for set B, and the 
stripe geometry, a trapezoidal cross-section of the stripe 
also shown in Fig. [4j The largest relaxation is observed 
in the corners of the stripes. 

The slopes of the sides in set B are few degrees larger 
than angles typically occurring when dry etching is used 
during the patterning process. Note that the X-ray 
diffraction reveals only the regions with regular lattice 
structure whereas the dry etching can leave a thin non- 
uniform amorphous coating on the stripes which leads to 
the unexpected non-rectilinear shape of the stripe cross- 
section resulting from the fitting. 

In the next step, we use our modelling of the lattice re- 
laxation also for stripes of set A where the X-ray diffrac- 
tion would be less accurate due to the small volume of 
the epilayer, however, the relaxation mechanism should 
be of the same nature as in set B. Fig. [5] shows the spa- 
tial dependence of the function p(x' , z') for two different 
geometries relevant to samples in set A. The shape of 
the stripe cross-section cannot be determined from the 
SEM image of Fig. [I] with nanometer accuracy. This 
uncertainty cannot be neglected in the undercut stripes 
A[no]. Therefore, more geometries (slopes of the sides) 
were simulated and one representative example is given 
in the upper panel of Fig. [5] On the other hand, the 
precise shape of the sides does not play such an impor- 
tant role in case of the overcut stripes Am^qi shown in the 
lower panel of Fig. [5] In all geometries, the local induced 
strain is stronger closer to the edges of the stripes. 

The comparison of the macroscopic simulations and X- 




FIG. 4: (Color online) Finite element calculation of the lat- 
tice relaxation, p(x',z'), on the cross-section perpendicular 
to the slightly overcut stripes Br.^ (upper panel) and B[ 110 ] 
(lower panel). The cross-section of one stripe and the under- 
lying substrate is plotted. The relaxation p = 1 and p = 
corresponds to a full relaxation of the lattice and to a lattice 
under a compressive strain of the unpatterned layer, respec- 
tively. Both stripes are close to full relaxation. 
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FIG. 5: (Color online) Finite element calculation of the lat- 
tice relaxation, p(x', z'), on the cross-section perpendicular to 
the undercut stripes A[ 110 ] (upper panel) and overcut stripes 
A [iioi (lower panel). The cross-section of one stripe and the 
underlying substrate is plotted. The relaxation p = 1 and 
p — corresponds to full relaxation of the lattice and to a 
lattice under a compressive strain of the unpatterned layer, 
respectively. All stripes show weaker net relaxation than the 
stripes in set B. 



ray diffraction measurements are done on the level of the 
full spatial distribution of the relaxation p(x',z'). The 
magnetic characteristics, considered in this work in the 
single domain approximation, are analyzed based on the 
net lattice relaxation. Here we take advantage of the di- 
rect proportionality of the magnetocrystalline anisotropy 
to the corresponding strai n 18 ' 19 ' and calculate the mean 
anisotropy from the spatial average of p(x' ,z') over the 
stripe cross-section. We will denote this average quantity 
by p in the rest of the paper. 

The last step in obtaining the input parameters for the 
microscopic modelling is writing the net in-plane compo- 



nents of the total strain tensor introduced in Eq. ([I]) in 
the coordinate system fixed to the main crystallographic 
axes introduced in Fig. [lj 
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e ( 1 - - - - sin 2u ) , 
e ( 1 - - + - sin 2uj ) . 



(2) 



e xy = e -cos2uj, 

where the angle uj is measured from the [110] axis and 
the angle uj — 45° describes the rotation of the relax- 
ation direction (the dashed coordinates) with respect to 
the crystalline coordinate system. Note that the above 
strain components coincide with those in Eq. (IT]) when 
uj = 45°, i.e., the relaxation direction is aligned with the 
[100] axis. We emphasize that the average relaxation p 
depends on uj. We rotate the elasticity matrix describing 
the cubic crystal when simulating the lattice relaxation 
along different directions. 

The strain components e XXl e yyi and e xy for the stripes 
in set A are obtained from th e mac roscopic simulations 
and considering eo ~ — 0.22%J- 16 * 2Q I Table y] summarizes 





eo[%] 


P 


A - 

^[110] 


-0.22 ±0.03 


0.184 ±0.005 


A[no] 


-0.22 ±0.03 


0.24 ±0.05 


B [iToi 


-0.38 ±0.03 


0.79 ±0.01 


B[oio] 


-0.38 ±0.03 


0.99 ±0.01 



TABLE I: The lattice mismatch eo and the lattice relaxation 
p for the patterned samples as entering the microscopic cal- 
culations in Sec. ^ The value of eo in set B is determined 
from the X-ray diffraction experiment, whereas eo in set A is 
inferred from the partial Mn concentrations using the analysis 
of Refs. [20TT6] . 

the parameters determined in this section. 



IV. EXPERIMENTAL MAGNETIC 
ANISOTROPIES 

In-plane magnetic anisotropics in thin (Ga,Mn)As 
films are often analyzed using the lowest order decom- 
position of the free energy profile into separate terms of 
distinct symmetry! 21 " 23 ! In this study, we follow this track 
by adopting the following phenomenological formula: 

F(i/>) = -^ sin 2 2ip ± K u sin 2 ip - K Q sin 2 (^ - ft). (3) 

The cubic symmetry of the underlying zinc-blende struc- 
ture is described by the first term with minima along the 
[100] and [010] directions in case of K c > 0. The sec- 
ond term quantified by the coefficient K u describes the 
so called "intrinsic" uniaxial anisotropy along the crystal 
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FIG. 6: (Color online) Remanent magnetization along the 
main crystallographic directions for sample Ao (25 nm thick 
unpatterned epilayer). 



diagonals present in the unpatterned (Ga,Mn)As epilay- 
ers. The last term quantified by Kq describes the uniax- 
ial anisotropy with an extremum at an angle Q induced 
by the relaxation of the lattice mismatch of the doped 
epilayer and the substrate. The angle Q is in general not 
equal to the angle of the corresponding lattice relaxation 
cjl^Both angles are measured from the [110] axis. 



A. Remanent magnetization 

Remanent magnetization along the main crystallo- 
graphic directions was measured by SQUID for both sets 
of samples. The obtained values include the magneto- 
crystalline anisotropics described in the previous para- 
graph as well as the shape anisotropy which always 
prefers the magnetization alignment with the longest side 
of a rectangular prism such as the stripesP^ 

Fig. [6] shows that in the control sample Aq the intrinsic 
uniaxial anisotropy dominates over the cubic anisotropy 
on a large temperature range and the easy axis along 
the [110] diagonal. The ratio of the remanent magneti- 
zation projections to the [110] and [100] directions below 
60 K reveals that the system is almost purely uniaxial. 
The behavior of the anisotropy components at T > 60 K 
cannot be described within the single domain approxima- 
tion. However, the anisotropics of unpatterned samples 
are relevant to our microscopic analysis of measurements 
in the microbars only at the lowest temperatures where 
we extract intrinsic anisotropy coefficients and deduce 
the material parameters as detailed in Sec. [V| 

Fig. [7| shows that the patterning of the sample Am^ 
strengthens the uniaxial anisotropy present in the par- 
ent wafer. The [110] diagonal becomes the easiest of the 
investigated directions at all temperatures and the [110] 
diagonal becomes the hardest axis at all temperatures 
below Tc- 

Fig. [8| shows that in the sample A[ 110 ], the two diago- 
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FIG. 7: (Color online) Remanent magnetization along the 
main crystallographic directions for sample Ar^i (750 nm 

wide stripes along the [110] direction). 
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FIG. 8: (Color online) Remanent magnetization along the 
main crystallographic directions for sample A[ 110 ] (750nm 
wide stripes along the [110] direction). 



nals switch roles and in analogy with the previous case 
the easy axis prefers alignment close to the stripe direc- 
tion, which is the hard axis over most of the temperature 
range in the parent wafer. This means that a rotation 
of the easy axis by as much as 90° is achieved by the 
post-growth patterning. Note that the difference of the 
projection of the remanent magnetization to the [100] 
and [010] directions in the two patterned samples is due 
to a 5° misalignment between the stripes and the crystal 
diagonals introduced during the fabrication. 

The samples in set B posses stronger cubic anisotropy. 
Fig. [9] shows that in the control sample Bq the intrinsic 
uniaxial anisotropy dominates over the cubic anisotropy 
only at temperatures above 20 K and the [110] diagonal is 
easier than the [110] diagonal at all temperatures below 

Tc- 

Fig. 10 shows a strengthening of the uniaxial 



anisotropy along the stripe direction in the sample 
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FIG. 9: (Color online) Remanent magnetization along the 
main crystallographic directions for sample Bo (200 nm thick 
unpatterned epilayer). 
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FIG. 10: (Color online) Remanent magnetization along the 
main crystallographic directions for sample Brjj i (1 fim wide 

stripes along the [110] direction). 



anisotropy at the lowest temperatures. The transition 
from cubic to uniaxial anisotropy occurs at a lower tem- 
perature than in the control sample. The [110] direction 
is hardened. The main crystal axes [100] and [010] re- 
main equal due to the more accurate alignment of the 
stripes with the crystal diagonal. 
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[110] 



although not large enough to overcome the cubic 



Fig. 11 shows a differentiation of the [100] and [010] 
projections in the sample Bpio]- The uniaxial anisotropy 
along the stripe direction now dominates at all tempera- 
tures. The intrinsic anisotropy differentiating the diago- 
nal directions is less pronounced than in case of Bo as it 
has to compete also with the induced uniaxial anisotropy. 

We can conclude that the universal effect seen in all 
patterned (Ga,Mn)As/GaAs samples is the preference 
of the easy axis to align parallel to the stripe which 
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FIG. 11: (Color online) Remanent magnetization along the 
main crystallographic directions for sample Bpio] (1 M m wide 
stripes along the [010] direction). 



is the direction in which the growth induced compres- 
sive strain cannot relax, i.e., the direction of the rel- 
ative lattice contraction in (Ga,Mn)As. This is remi- 
niscent of the magnetocrystalline anisotropy of unpat- 
terned (Ga,Mn)As epilayers which typically yields easy- 
axis oriented also along the direction of contraction, i.e., 
in-plane for compressively strained (Ga,Mn)As epilayers 
and out-of-plane for (Ga,Mn)As films grown under ten- 
sile strain. 19 We point out that the measured magnitudes 
of magnetic anisotropics in the microbars are an order 
of magnitude larger than the shape anisotropy contribu- 
tion for given concentration of magnetic moments and 
thickness to width ratio. The microfabrication effects 
in the (Ga,Mn)As stripes are therefore primarily due 
to the spin-orbit coupling induced magnetocrystalline 
anisotropy. 



B. Anisotropy coefficients 

After investigating the reorientations of the easy axis 
we focus on the magnitude of the individual anisotropy 
components. We measure the hysteresis loops using the 
SQUID magnetometry and fit the results to the following 
equation: 
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(4) 



-M s H n sin 2 (/0 - £2) - M s Hcos(ip - <j> H ) 



where K{ = noMsHi were introduced in Eq. (J3]) , Ms is 

the saturation magnetization, H is the external magnetic 
field applied at the angle 4>h, and the last term is the Zee- 
man energy. All angles in Eq. Q are measured from the 
[110] axis. In case of a general alignment of the induced 
uniaxial strain, the angle Vt of the corresponding uniaxial 
anisotropy is an independent fitting parameter. However, 



in case of the main crystallographic axes and their small 
surrounding we can set O = u;, i.e., the anisotropy term 
is aligned with the corresponding uniaxial strainP^ An 
overview of the resulting angles Q for the different align- 
ments of stripes in sets A and B is given in Table ITTl 





K c [kJ/m 3 ] 


K u [kJ/m 3 ] 


K n [kJ/m 3 ] 


Q [deg] 


A 
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0.404 
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A[no] 


0.412 


0.404 
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Bo 


2.213 


0.381 


0.0 




B [iTo] 


2.213 


0.381 


0.935 


90 


B[oio] 


2.213 


0.381 


0.696 


45 



TABLE II: The anisotropy coefficients obtained by fitting the 
hysteresis loops at T = 2 K to Eq. (El and the angular shift 
of the anisotropy term induced by the lattice relaxation as 
introduced in Eq. (13b - Note that the lattice relaxes perpen- 
dicular to the stripe direction. The error of the anisotropy 
coefficients is approximately 10 — 20%, approaching the up- 
per limit in case of the thick inhomogeneous samples in set B. 

When determining the anisotropy coefficients in the 
stripes we use the assumption of linear superposition of 
the anisotropics present in the unpatterned samples with 
the anisotropics induced by the patterning and lattice 
relaxation: the coefficients K c and K u are obtained first 
in the control samples and kept fixed when fitting the 
stripes fabricated from the same epilayer. The assump- 
tion is justified on the qualitative level by the rema- 
nent magnetization measurement discussed in the pre- 
vious subsection which revealed the persistence of the 
bulk anisotropics in all patterned samples. Its valid- 
ity has been corroborated also by studies of epilayers 
subject to post-growth piezo straining 3 and lithographic 
patterningP We emphasize that our approach is appro- 
priate only when the lattice relaxation direction is very 
close to the main crystallographic axes or when the angle 
Q is also treated as a fitting paramaterP^ 

Another assumption concerns the magnetization reori- 
entation mechanism determining the shape of the hys- 
teresis loops. In case of a dominant uniaxial anisotropy 
we fit the hysteresis loops obtained for external fields ap- 
plied along the hard axis. In case of a dominant cubic 
anisotropy there is no completely hard direction. We 
nevertheless still consider a single domain model in the 
fitting. 

Anisotropy coefficients for all six samples at the low- 
est temperature are summarized in Table [TTJ Recall that 
these energies include also the contribution of the shape 
anisotropy which amounts to ~ 0.1 kJ/m 3 in samples 
A[iToi an d A[no] and ~ 0.3 kJ/m 3 in the samples Bm^qi 
and B[oio] with the higher thickness to width ratio. Note 
that the smaller coefficient K^ leads to the formation of 



a strongly uniaxial system as shown in Figs. 11, whereas 
the larger coefficient Kqq cannot overcome the cubic 
anisotropy component, at least at low temperatures as 




FIG. 12: (Color online) Anisotropy coefficients as functions of 
temperature obtained by fitting the hysteresis loops to Eq. (|3| 
for the three samples of set A. The uniaxial coefficients Kn 
(denoted by Ka for set A) due to the growth strain relaxation 
in the patterned samples dominate the total anisotropy. 
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FIG. 13: (Color online) Correspondence of the hole density, 
p, and the intrinsic shear strain, e^n, to the effective Mn con- 
centration, x, based on the agreement of the calculated K c 
and K u with the measured values. Samples A and B at zero 
temperature. 



shown in Fig. 10 It is because in case of sample Bjqio]? 
the induced anisotropy is added along the [010] axis 
which was already the easy (together with [100]) direc- 
tion in the unpatterned epilayer. 

For the thinner and more homogeneous epilayers in 
set A we were able to extract the temperature depen- 
dence of the anisotropy coefficients from the hysteresis 
loops up to T = 60 K as shown in Fig. [12] The uni- 
axial coefficients due to lattice relaxation dominate the 
anisotropy at all temperatures. At low temperatures the 
relative size of the induced anisotropics corresponds well 
to the simulated relaxations p: Sample Ar^oi with over- 
cut sides (weaker relaxation) shows smaller anisotropy 
than sample A[ 110 ] with undercut sides (stronger relax- 
ation). The cubic anisotropy remains positive for all 
studied temperatures T < 60 K which is in good agree- 
ment with the remanent magnetization data shown in 
Fig. [6j We do not discuss measurements above 60 K for 
which, as mentioned above, the single domain model is 
not applicable. 



V. COMPARISON WITH THEORY 



In this section we build on macroscopic calculations 



of the lattice relaxation presented in Sec. Ill, perform 
the microscopic calculations of the magnetic anisotropy 
energy, and analyze its correspondence with the experi- 
mental results on the level of anisotropy coefficients. We 
extract the coefficients by fitting the calculated total en- 
ergies to Eq. pi) for different magnetization directions. 

The comparison of the experimental and theoretical 
results involves a number of material parameters. The 
most important inputs of the microscopic calculations 
are the concentration of the ferromagnetically ordered 
Mn local moments (x) and the hole density (p). Unfor- 
tunately, these two parameters cannot be accurately con- 



trolled during the growth or determined post growthPSl 
The measured saturation magnetization, the conductiv- 
ity, and the Curie temperature of the control samples 
provide only estimates of these input parameters with 
limited accuracy. 

Another independent input parameter of the micro- 
scopic simulations is the "intrinsic" shear strain which 
has been used successfully to modeP^^the in-plane uni- 
axial anisotropy in the unpatterned samples. We re- 
call that such modelling for small strains (the typical 
values 19 are e^ n ~ 10 -4 ) complies well with the assump- 
tion that the "intrinsic" uniaxial anisotropy superposes 
linearly with anisotropy components induced by the lat- 
tice relaxation, as mentioned in the previous section. The 
intrinsic shear strain is added to the off-diagonal ele- 
ment of the total strain tensor written in Eq. (|2| giving: 
e xy = e in + e f cos 2a;. 

Fig. 13 shows the combinations of x, p, and e^ n for 
which the calculated intrinsic K u and K c of the con- 
trol samples Ao and Bo agree with the measured values 
at zero temperature. By this we limit the intervals of 
x, p, and ei n values considered in the modeling of the 
temperature dependent anisotropy coefficients in all mea- 
sured samples. Note, that we have also imposed an upper 
bound to x given by the nominal Mn doping in the par- 
ticular material and a bound to p ensuring a maximum of 
one hole per Mn ion and in-plane easy axis (axes). This 
method allows for predicting the induced anisotropy co- 
efficients in the microbars without any adjustable param- 
eters in the microscopic model. 



A. Low temperatures 



Using parameter combinations shown in Fig. 13 we cal- 
culate the induced uniaxial anisotropics in the microbars 
at zero temperature. The left and right vertical axis of 
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FIG. 14: (Color online) Calculated anisotropy coefficients due 
to the lattice relaxation in the patterned samples A and B at 
zero temperature for fixed combinations of x and p shown in 
Fig. pi 
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FIG. 15: (Color online) Calculated cubic a uniaxial intrinsic 
anisotropy coefficients present in all samples A as functions 
of temperature for fixed combinations of x, p, and ei n shown 
in Fig. fT3l 



Fig. [14] shows the extracted anisotropy coefficients for 
stripes in sets A and B, respectively. The combinations 
of x, p, and 6i n are indexed only by x for simplicity. The 
plotted values can be compared to the measured coeffi- 
cients summarized in Table [TTJ The relations K^ < K$ 
and Kgo > K45 hold both in theory and in experiment. 
We observe a semi-quantitative agreement in samples 
AfiToi an d A[no] where the measured values are roughly 
a factor of 2 larger than the calculated ones. The ra- 
tio of the calculated coefficients Kq for samples Am^ 
and A[no], K95/K5, is in excellent agreement with ex- 
periment (the difference is only 4%). These agreements 
justify the interpretation of the measured effects in the 
microbars based on the strain-relaxation controlled mag- 
netocrystalline anisotropy. Note that they also support 
the assumption of the linear superposition of individual 
uniaxial anisotropics terms used in our analysis. 



Fig. 14 shows also extracted anisotropy coefficients for 
samples B^^ and B[ io]- I n this case, the calculated 
ratio of coefficients extracted for the two stripe align- 
ments, If 90 /if 45, is approximately 20% larger then the 
corresponding experimental ratio, i.e., still in a very good 
agreement. We note, however, that the absolute values 
of the measured coefficients are about a factor of 10 lower 
than the calculated ones. A possible source of this dis- 
crepancy is the large value of the experimentally inferred 
K c due to inaccurate subtraction of the paramagnetic 
and diamagnetic backgrounds from the measured hys- 
teresis curves. In general, we also expect that the theo- 
retical modelling is less reliable in the thicker, as-grown 
samples B due to stronger disorder and inhomogeneities 
in the material. 

As a consequence of the almost complete relaxation 
of the lattice mismatch in the thicker samples the cal- 
culated anisotropy coefficients are larger than the cubic 
coefficient at all studied temperatures which is not in 
agreement with the measured coefficients in set B at low 
temperature (see Table [ll|) . 



B. Temperature dependence 

We now select six representative combinations of x, 
p, and ei n from the relevant interval shown in Fig. [l3j 
calculate the temperature dependence of all anisotropy 
coefficients for each set of parameters, and discuss the 
comparison with the measured anisotropics. We recall 
that in our mean-field modeling at finite temperatures 
the calculated Tc is uniquely determined by x and p. 
Note that for the entire interval of relevant x and p deter- 
mined from the low-temperature analysis in the previous 
section, we obtain Curie temperatures which are in agree- 
ment with the experimental values in materials A and B 
within a factor of 2. This provides an additional support 
for the overall consistency of our microscopic theoretical 
analysis of the measured data. 

Fig. [15] shows the calculated intrinsic anisotropy coef- 
ficients K c and K u of samples in set A for three fixed 
parameter combinations. At zero temperature the values 



coincide with data in Fig. [I3j The cubic anisotropy com- 
ponent is stronger than the intrinsic uniaxial component 
at lowest temperatures but it quickly becomes weaker 
as temperature is increased for all parameter combina- 
tions. This temperature dependence is in agreement with 
the experimental anisotropics measured below 60 K, as 
shown in Fig. 12 The comparison cannot be extended 



to higher temperatures because, as explained above, the 
experimental behavior at these temperatures is not cap- 
tured by the single domain model. 

Fig. [16] shows the calculated anisotropy coefficients K^ 
of samples An^ and A[ 110 ] again for the three fixed pa- 
rameter combinations. The calculated anisotropy com- 
ponents decrease monotonously with increasing temper- 
ature in agreement with the measured dependencies pre- 
sented in Fig. [T2l The comparison provides additional 



support for the interpretation of the experimental data, 
suggested already by the analysis at low-temperature, 
which is based on the strain relaxation induced magne- 



10 



0.7 

0.6 

0.5 

"S 0.4 

< 0.3 

0.2 

0.1 

0, 





1 ' 1 ' 1 ' 1 


1 


I 








A Mtm , x = 3% 


- 








^•^ 




[110]' 






S 





A [-110]' X = 3% 




=-= 







A [U0] ,x = 4% 










A r 1im , x = 4% 






\ S N 




[-110]' 


- 


- 


— \ V N 





A [U0] ,x = 5% 


- 


- 






A [ . U0] ,x = 5% 


- 










A s A 










A w 


\ 




\ 


'•*.s 




V N ^ 


■•>.• 




i 1 i 1 i > I i 1 rv. 


1 1 V | 



20 



40 



60 



T[K] 



100 



120 



140 




B, 110] ,x = 5% 
B [010] ,x = 5% 



■■ B [0 ior x = 6% 

B, 110] ,x = 7% 

■■ B [oior x = 7% 



100 120 140 160 180 
T[K] 



FIG. 16: (Color online) Calculated anisotropy coefficients due 
to the lattice relaxation in the patterned samples A as func- 
tions of temperature for fixed combinations of x and p shown 
in Fig. [13] and for the induced strain given in Fig. ^4] 



FIG. 18: (Color online) Calculated anisotropy coefficients due 
to the lattice relaxation in the patterned samples B as func- 
tions of temperature for fixed combinations of x and p shown 
in Fig. [13] and for the induced strain given in Fig. ^4] 
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FIG. 17: (Color online) Calculated cubic a uniaxial intrinsic 
anisotropy coefficients present in all samples B as functions 
of temperature for fixed combinations of x, p, and e^ n shown 
in Fig. O 



tocrystalline ansisotropy effects. 

Fig. [Tt] shows the calculated intrinsic anisotropy coef- 
ficients K c and K u of samples in set B again for three 
fixed parameter combinations. At zero temperature the 



values coincide with data in Fig. [13] The calculated cubic 
anisotropy dominates over the uniaxial anisotropy at low 
temperatures in agreement with the experiment. The 
cross-over in the theory curves to the dominant uniax- 
ial anisotropy occurs at higher temperatures than Tc/3 
observed in experiment (see Fig. |9|; at the upper part 
of the relevant interval of Mn concentrations the theo- 
retical crossover occurs at Tq/2. We again attribute this 
quantitative discrepancy to inhomogeneities and stronger 
disorder in the thick as- grown material B. 



Fig. 18 shows the anisotropy coefficients Kq of sam- 
ples Bmjqi and B[ io] for the same fixed parameter com- 



binations as in Fig. 17 Again, we observe very simi- 
lar dependence of the uniaxial anisotropy coefficients on 
temperature as in experiment. The monotonous decrease 
of the coefficients with growing temperature is in agree- 



ment with the measured remanent magnetization data 
in Figs. [TT] and [To] Both induced anisotropy coefficients 
are predicted to be larger than the cubic coefficient at all 
studied temperatures. This complies with the measured 
remanence data of sample Bpio]- Sample Bm^qi shows 
agreement above 20 K. Its behavior at temperatures be- 
low 20 K, is not captured by the theory data as we have 
already discussed in the previous subsection. 



VI. SUMMARY 

We have performed a detailed experimental and theo- 
retical analysis of magnetic anisotropics induced in litho- 
graphically patterned (Ga,Mn)As/GaAs microbar ar- 
rays. Structural properties of the microbars have been 
studied by X-ray spectroscopy showing strong strain re- 
laxation transverse to the bar axis. The relaxation in- 
duced lattice distortion in stripes with thickness to width 
ratio as small as ~ 0.1 induces additional uniaxial mag- 
netic anisotropy components which dominate the mag- 
netic anisotropy of the unpatterned (Ga,Mn) As epilayer, 
as revealed by SQUID magnetization measurements. The 
easy axis can be rotated by the micropatterning by 90° 
at all temperatures below the Curie temperature. 

We have carried out systematic macroscopic and mi- 
croscopic modeling of the structural and magnetic char- 
acteristics of the microbars and analyzed in detail the ex- 
perimental results. The agreement of the measured and 
simulated X-ray diffraction maps shows that the applied 
elastic theory model is quantitatively reliable in predict- 
ing the local lattice relaxation in patterned epilayers with 
the growth induced strain. The overall good agreement 
of the microscopically calculated and measured magnetic 
anisotropics conclusively demonstrate that the pattern- 
ing induced anisotropics are of the magnetocrystalline, 
spin-orbit coupling origin. 
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